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Abstract — The coupling between negative and positive index 
medium waveguides is investigated theoretically in this paper. A 
coupled mode theory is developed for such a waveguide system 
and its validity is verified. Interesting phenomena in the coupled 
waveguides are demonstrated, which occur in the case when the 
negative index medium waveguide in isolation guides its mode 
backward. A new type of coupled mode solution that varies 
exponentially with the coupling length is found in the special case 
when the propagation constants of two individual waveguides are 
nearly the same. A coupler operating in this case is insensitive to 
the coupling length, and its coupling efficiency can reach 100% 
as long as the coupling length is long enough. However, when 
the propagation constants of the two individual waveguides differ 
greatly, the coupled mode solution is still a periodic function of 
the coupling length, but the coupled power is output backward. 
In addition, the modes in the composite waveguide system are also 
studied using the coupled mode theory, and their fundamental 
properties are revealed. 

Index Terms — Negative index media, waveguides, coupled 
mode theory, backward waves. 



I. Introduction 

Negative index media (NIM) that have simultaneously neg- 
ative permittivity and permeability, have attracted intensive 
interest recently. These media exhibit several extraordinary 
effects such as negative refraction, backward waves, and 
evanescent wave amplification [l]-[3], seeming to challenge 
several concepts well established for familiar positve index 
media (PIM) in electromagnetism and optics. NIM materials 
are not generally found in nature, and thus need to be artifi- 
cially constructed. Though the concept of NIM was already 
proposed by Vesalogo in 1968 [1], it was not until 2001 
that the first experimental demonstration of negative index 
behavior was accomplished by Shelby et al., with a material 
made by a two-dimensional array of repeated unit cells of 
copper strips and split-ring resonators [4]. So far, a variety of 
NIM structures have been proposed [4]-[8] in the microwave 
regime. More recently, NIM structures operating in the infrared 
and visible frequency range have also been reported [9]— [13], 
e.g., a NIM formed by an array of pairs of parallel gold rods 
was demonstrated to have a negative refraction index at a 
wavelength of 1.5 /im [10]. 
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While the study of various NIM structures has been a 
subject of growing interest during the past few years, several 
research works have involved waveguides with NIM compo- 
nents, and interesting properties of guidance were reported 
[14]— [16]. For a NIM slab waveguide, the portion of a guided 
mode inside the slab has a Poynting vector contradirectional 
with the direction of the phase velocity of the mode, while the 
portion of the guided mode outside the slab has a Poynting 
vector parallel to the phase velocity. The total energy flow 
of the guided mode may be codirectional or contradirectional 
with the phase flow. No fundamental mode exists in a NIM 
slab, but it may still support a single-mode propagation under 
certain conditions. Moreover, the interaction between PIM 
and NIM waveguides was investigated in [17], and the phe- 
nomenon of anti-directional coupling was demonstrated, which 
offers a new possibility in the design of future devices and 
components. The bound modes of a NIM and PIM composite 
waveguide were also studied in [18], [19]. 

The coupling of waves in planar dielectric waveguides 
provides a means of reflectionless signal transfer from one 
waveguide to another, thus plays an important role in inte- 
grated optics. Compared to the conventional waveguide sys- 
tem [20]-[22], the coupled NIM and PIM waveguide system 
receives less attention. In this paper we will study carefully 
the coupling between planar NIM and PIM waveguides. For 
this purpose, a coupled mode theory is developed for such 
a waveguide system. Our analysis will show that there exist 
various possibilities of coupled mode solution for the NIM 
and PIM waveguide system, and the interaction of the guided 
modes of individual waveguides may even be equivalent to the 
interference of evanescent modes in the composite waveguide 
for particular cases. In Section II, the coupled mode theory for 
parallel NIM and PIM waveguides is formulated, and various 
coupled mode solutions are discussed in Section III. In Section 
IV, the guidance properties of a NIM and PIM composite 
waveguide are analyzed. Also, the validity of the coupled 
mode theory is examined in this section. Section V concludes 
the paper. 

II. The Coupled Mode Equations 

Consider two parallel planar waveguides as illustrated in 
Fig. 1. The upper waveguide layer of width a is a PIM with 
relative permittivity e a > and permeability /i a > 0, and 
the lower one of width b is a NIM with relative permittivity 
£b < and permeability < 0. The core layers are 
separated by a distance of d and they are surrounded by 
free space. We assume that the individual PIM waveguide 
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Fig. 1. Schematic of the coupled waveguide structure. 



supports only one guided mode with electric and magnetic 
fields {ei(x)exp(jPiz),hi(x)exp(jP 1 z)} (/3\ is the propa- 
gation constant of the guided mode), whereas the individual 
NIM waveguide may allow multi-mode propagation. The 
modal fields of the individual NIM waveguide are denoted by 
{e m (x) exp(jf3 m z), h m (x) exp(j/3 m z)} (J3 m are the propaga- 
tion constants of the modes), where m = 2, 3, • • -,N. We 
consider an electromagnetic field {~E(x, z),H(x, z)} which 
satisfies Maxwell's equations plus the boundary conditions of 
the entire structure. Its transverse field can be written as 

N 



(x)] , (1) 

m— 1 
N 

Hi (a;, z) = ^ u m(z) \hmt(x) + Sh m (x)} , (2) 



where 5e m and 5h m (m = 1, 2, • • •, N) are the corrections 
of the modal field in the other guide due to the induced 
polarization. The correction fields Se m are introduced as in 
Ref. [22], and <5h TO are treated in a similar manner. The 
longitudinal components of the vectorial field are expressed 
in the form 



N 



E 2 (x, z) = V* u m (z) £m } X } e mz (x), (3) 

z — ' E[X) 



m—1 

N 



n z ( X , Z)=J2 "mW^hm.W, . (4) 



where e\ — e'"', 
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2. {e,^}, {e(°),^ )}, and {e^ 



evi, ix m = /i( b ) for m > 
/i^} represent the profiles 
of the entire structure, and of the individual PIM and NIM 
waveguides, respectively. 

We utilize a variational principle for the coupled system 
under consideration following Haus et al. [21], and obtain the 
differential equations in the form 



P—U = -jPBU-jKU, 
dz 



(5) 



K 
with 



and 

u± ran 
AP 

SK mn 
AK 



{SP mn + SP* m )+AP mn , (6) 
- SK mn + ([3 m - p n )SP mn + AK mn , (7) 



[e* mt x h nt + e nt x h* mt ] ■ zdx, (8) 

J[Ae n e* n ■ e„ - {Ae m Ae n /e)e zn ■ e* zm 

A/x n h^ • h„ - (AfimAfin/ fi)h. zn ■ h* zm ]dx, 

(9) 

J [5e n x h* nT + e* mT x Sh n ] ■ zdx, (10) 

J [Se* n x Sh n + Se n x <5h^J • zdx, (11) 

+ Ae„i5e* n • e„ + A/in^h^ • h„ 
A^i m h^ • <5h„]da;, (12) 



[eSe* m ■ 5e n + ^Sh* m ■ £h n + (V t x <5h„) ■ Je* 



where C/ = colfui, «2, • ■ •, mat], and B = diag[/?i, $2, ■ ■ ■, /?jv]- 



- (V t x 5e„) • Sh* m ]dx, (13) 

where Ae m = e — e m and A/i m = /i — /i m . Note that AP„ m = 
(5^mn = Aif„ m = for m^n. 

The above coupled-mode equations are similar to those 
derived by Haus et al. in [23], but our formulation considers 
a more general case, in which the permeability of the coupled 
waveguides is allowed to vary with the transverse coordinate 
(x). Furthermore, our formulation includes the factor Eije 
(/J-i/n) in the expression for E z (H z ). We let H = PB+K = 
H + SH, here H = PB + K, then Eq. (5) can be written 
as PdU/dz = —jHU. In the case if the waveguide system 
is lossless, H can be proved to be Hermitian in the same 
manner as in [10]. As the matrix SH has off-diagonal elements 
SH mn = (3 m 5P mn + f3 n SP* m , thus the matrix H is also 
Hermitian for the lossless case. So the power conservation 
holds for our formulation [24]. 

III. Coupled Mode Solutions 

Let us consider the wave propagation and coupling along 
two lossless planar waveguides over a finite length L (see 
Fig. 1). For simplicity, we assume that the individual NIM 
waveguide also supports a singled mode (i.e., N = 2). 
The single-mode operation for a certain polarization in the 
NIM waveguide can be achieved through suppressing the 
appearance of surface modes. We consider an initial power 
to be injected into the PIM waveguide at z = 0. The excited 
mode in the NIM waveguide will have the same direction of 
phase flow as the guided mode in the PIM waveguide, i.e., 
02 > 0. As the energy flows of the mode in the core and 
cladding of the NIM waveguide are in the opposite directions, 
the total energy flow may be codirectional or contradirectional 
with the phase velocity of the mode, i.e., P22 > or P22 < 0. 
The coupled power in the NIM waveguide will be output at 
the waveguide end z = L if P22 > 0. Otherwise the coupled 
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power will be output at the other end z = 0, if P22 < 0. 
Correspondingly, the boundary condition is u 2 = at z = 
for P22 > 0, and it becomes u 2 — at z = L for P22 < 0. 
We rewrite Eq. (5) in the form 

d 



dz 



U = -jBU - jCU, 



(14) 1 



where C = P~ 1 K = P~ X H — B, with the superscript 
"— 1" denoting an inverse matrix. The coefficients C mm 
(to = 1, 2) in Eq. (14) are equivalent to a modifica- 
tion of the parameters (3 m , and our analysis shows that 
Cmm I ^ A™, in general. The coupling coefficients C12 and 
C21 are especially interesting. These coupling coefficients are 
given by C 12 = (P22P12 - P12H22) / {P11P22 - P12P21) 
and C21 - (Pxiflai - P2iH u ) / {P 1X P 22 - P12P22), where 
#11 = PnA + #11 and H22 = P22P2 + K 22 . Since 
IP12-P21I « I-P11P22I, \K 22 \ « LP 2 2/3 2 |, and P n « P u 
we find C12 ~ (P12 — Pi2/32)/Pn- Similarly, we obtain 
C21 w (#21 - P2i/3i)/P22- Note that P i2 = H 2 i and 
P12 = P21 in the lossless case. If /3i and (32 differ negligibly, 
i.e., (3i w /?2, then C12 and C21 have the same sign for 
P22 > 0, or they have opposite signs for P22 < 0, which 
never happens for coupled conventional waveguides. As an 
example, the coupling coefficients as a function of the NIM 
layer thickness b are illustrated in Fig. 2(a), and the value of 
(fii — Pi) is shown in Fig. 2(b). The guided wave is assumed 
to be TE-polarized, and the parameters of the two-waveguide 
system are as follows: E\ = 2.1, [i\ = 1, and a = 0.3035A, 
corresponding to the waveguide parameter V = 1 for the PIM 
waveguide; £2 = —1.8, and /12 = —1; d = 2a. As seen 
from Fig. 2(a), in the region of 0.36A < b < 0.55A where 
P22 > 0, the coupling coefficients are both positive; but they 
have opposite signs in the region of 1.12A < b < 1.63A where 
P22 < 0. The shaded areas in Fig. 2 indicate the regions 
where the single-mode propagation in the individual NIM 
waveguide is not available or only the surface mode exists. In 
what follows, we concentrate on the analysis of the coupling 
between PIM and NIM waveguides for TE polarization (results 
for TM polarization can be obtained from duality). We will 
consider two cases for waveguide system: (i) P22 < 0; (ii) 
P22 > 0. 

A. Case (i): P 22 < 

The solution satisfying u\ = ui(0) at z = and U2 = at 
z = L is shown to be 

ui(z) = A{^-j'Atanh[?A(z-i)]}e J02 , (15) 
u 2 (z) = jAC 2 itanh[V>(z-£)]e^ z (16) 

for the case if A 2 < — C12C21, or else it becomes 

m{z) = A{ipcos[ip(z-L)} - jAsin[^(z-£)]}e^ z , 

(17) 

u 2 (z) = AC 2 ism[^(z-L)]e^ z , (18) 

where V - y/\A 2 + C 12 C 2 i\, A = - 2 )/2, and 
<p=i$i + /3 2 )/2 with /3 m =P m + C mm (to = 1, 2). Since 
C12C21 < at least when Pi rj f3 2 (i.e., A w 0), the first 
type of solution, i.e., Eqs. (15) and (16), is possible to occur 
for the case with P22 < 0. Note that Eqs. (15) and (16) are a 




Fig. 2. (a) Coupling coefficients versus the NIM layer thickness b. (b) 
Propagation constant difference (@2 — Pi) versus the NIM layer thickness b. 
The shaded areas indicate the regions where the single-mode propagation in 
the individual NIM waveguide is not available. 



new type of the solution for coupled mode equations, which 
never happens for a conventional waveguide system. 

The coupling efficiency (or the fraction of the cou- 
pled power), which describes the net power transfer be- 
tween the two waveguides, is determined by r\ = 
|p22|u 2 (0) 2 /Pi«i(0) 2 - In the case of A 2 < -C 12 C 2U 
we have 77 = |P 22 |C| 1 /Pii[V' 2 coth 2 (-0i) + A 2 ], thus the 
output power by the NIM waveguide increases as the length 
L increases, and it approaches its maximum of r) max — 
IP22C21/P11C12I as L > I/VIC12C21I. Evidently, 
can reach nearly 100% when /3i « p 2 - In the other case if 
A 2 > — C12C21, where (3i and /?2 differ greatly, the coupling 
efficiency becomes 77 = |P 2 2 IC^/Pn [tp 2 cot 2 (ipL) + A 2 ]. In 
this case we find that 77 varies periodically with L, and it 
has a maximum of r) max = |p22C| 1 /PnA 2 | at the lengths 
of L = (n — l/2)ir/ijj, here n is an integer. Since A 2 > 
-C12C21, we have rj max < \P22C2\I P\\C\2\- For this case 
our numerical analysis indicates that r\ ma x is always less than 
1. To illustrate the typical coupling behaviors, Fig. 3 shows 



the fraction of the coupled power for TE modes as a function 
of L for the waveguide systems with b = 1.2947A and 1.228A. 
The other parameters of the waveguide systems are the same 
as in Fig. 2. The propagation constant of the individual NIM 
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waveguide has ft = ft for b = 1.2947A and ft = 0.9816ft 
for b = 1.228A, corresponding to two typical cases discussed 
above. Here, we should point out that a NIM is inherently 
lossy, i.e., both Cb and jU& have an imaginary part. Thus, the 
propagation constant (ft) of a guided mode in the individual 
NIM waveguide also has an imaginary part. In the case if the 
NIM is severely lossy, the new type of coupled mode solution 
will not occur even when ft is equal to the real part of ft. 
But such a NIM with considerable loss is not suitable for 
application in devices. By using low loss materials to structure 
NIM, it is possible to reduce the loss of a NIM to very low 
level at certain frequencies, which is the case of our interest. 




— P 2 = Pi 

P 2 = 0.981 6|3 1 
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Fig. 3. Fraction of coupled power 77 (%) versus the coupling length L. 



by 



B. Case (ii): P 22 > 

The solution with u\ = ui(0) and u 2 = at z = is given 

ui(z) = A[tpcos(ipz) - iAsiri(V>z)]e #z , (19) 
u 2 (z) = ACnwa^e***, (20) 



where <fi, A, and ip are defined above. In this case the coupling 
efficiency is determined by 77 = \P 22 \u 2 (L) 2 / PnUi(0) 2 , and 
it follows that r] = P 22 C 21 sin 2 (?/>L)/P 11 ?/> 2 , which is identi- 
cal to that for two parallel conventional waveguides. The max- 
imum of the coupling efficiency is r\ max — P 22 C 21 / Pnip 2 , 
occurring at the lengths of L = w(n — l/2)ip, where n is an 
integer, and it reaches nearly 100% when ft is equal to ft. 

IV. Guidance of the Composite Waveguide 

It is interesting to analyze the guidance characteristic of 
the composite two-core waveguide, as the coupling between 
guided modes of two individual waveguides may always be 
viewed as the beating of the modes of the composite structure. 
The propagation constants of the modes in the composite 
waveguide are easily derived from Eq. (14), and are given 

by 



a v / (/3;-/3 2 ) 2 +4C 12 C 21 

ftt = ^— ± g ' (21) 

where f3 m (m = 1, 2) are defined in the previous section. 
It is known that the exact dispersion relation for the planar 



composite waveguide can be derived analytically for both TE 
and TM polarizations. The electric field (E) for TE modes in 
the five regions of the waveguide structure (see Fig. 1) can be 
written as 

Ei = yA\ exp(— ax + i(3z), x > a (free space) 
E 2 = y [A 2 exp(ikpx) + B 2 exp(— ikpx)] exp(ift), 

< x < a (PIM) 
E3 = y [A3 exp(— ax) + B3 exp(aa;)] exp(zft), 

— d < x < (free space) 
E4 = y [A4 exp(ikNx) + B4 exp(— ik^x)] exp(ift), 

-(d+b) < x < -d(NIM) 
E 5 = yA 5 exp(ax + i(3z), x < —(d+b) (free space) 



w here a = y/(3 2 - k 2 , k P = y 1 e a [i a k 2 - f3 2 , k N = 
\J EbUbk 2 — (3 2 , and fco is the wave number in free space. The 
nonzero components of the magnetic field (H) can be obtained 
straightforwardly from E. The dispersion relation of TE modes 
is determined by imposing matching conditions on the parallel 
components of E and H at the interfaces x — — (b + d), —d, 
0, and a. Following this procedure yields 

[(/1 - V/i) -2cot(fc P a)][(/a - I//2) - 2 cot(k N b)]e ad 



= (/i+l//i)(/ 2 + l// 2 )e 



-<xd 



(22) 



with /1 = kp/(ctfj, a ) and f 2 = k]y/(afib)- This equation is 
the exact dispersion relation for TE modes in the composite 
waveguide. The exact dispersion relation for TM modes can 
be obtained by substitution of fj, «-> e in Eq. (22). However, 
it is very difficult to solve analytically the transcendental 
equation (22). In contrast, the formula (21) enables one to 
gain more insight into the characteristic of the waveguide 
system. In the following we focus on the analysis of TE 
modes in the composite waveguide (results for TM modes 
can be obtained from duality). We will show that Eq. (21) is a 
good approximation to Eq. (22). Without losing generality, the 
interacted modes in the individual PIM and NIM waveguides 
are assumed to have positive propagation constants, thus the 
real part of the propagation constant of each related mode in 
the composite waveguide is also positive. All results calculated 
from Eq. (21) will be compared with the exact results obtained 
from Eq. (22). This provides an approach to examine the 
validity of the coupled mode theory. 

We first consider the waveguide system with P 22 < 0. In the 
special case with ft « ft, where C\ 2 C 2 \ < 0, one sees from 
Eq. (21) that the propagation constants have an imaginary part 
and be conjugate for any separation distances d. Note that both 
the NIM and PIM in the waveguide system are assumed to be 
lossless here. So these modes in the composite waveguide are 
evanescent waves. The appearance of the evanescent modes 
is only a consequence of the special coupling between the 
PIM and NIM waveguides, which corresponds to the new 
type of coupled mode solution. In the case if ft and ft 
differ substantially, the propagation constants may also have 
an imaginary part for small d, namely, when the coupling 
between the individual waveguides is strong. However, as d 
is so increased that \C\ 2 C 2 \\ < (ft — ft) 2 , the propagation 
constants will become real. In this case both modes in the 
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composite waveguide are propagating waves. In the case if j3-y 
and 02 differ greatly, the propagation constants are always 
real for any d. To illustrate these behaviors, Fig. 4 shows 
the propagation constants for the composite waveguides with 
different NIM layer thicknesses b = 1.2947A, 1.228 A and 
1.1559A. The other parameters are the same as in Fig. 2. 
The propagation constants of the individual NIM waveguide 
are /3 2 = /3i, 0.9816/?i, and 0.9531/3i for the three cases, 
respectively. In Fig. 4, the exact results (labelled by "Exact") 
for the three cases are also included for comparison with the 
results (labelled by "CMT") obtained from Eq. (21), and the 
excellent agreement is observed. 
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Fig. 4. Propagation constants (3± versus the distance d between the core 
layers for the composite waveguides with various NIM layer thicknesses. The 
NIM waveguide in isolation guides a single mode backward (P22 < 0) for 
all cases, (a) Real parts of f3± ; (b) imaginary parts of f3± . The lines represent 
the results obtained from coupled mode theory, and the marks represent the 
results calculated from Eq. (22). 



Next, we consider the waveguide system with P22 > 0. In 
this case, our numerical analysis indicates that C12C21 > 0, 
and the propagation constants of the modes of the composite 
waveguide are always real, as illustrated in Fig. 5, where the 
thicknesses of the NIM layers are b = 0.3832A, 0.3895A and 
0.4007A, corresponding to the propagation constants of the 
guided modes in the individual NIM waveguides /3 2 = Pi, 
0.9816/?i, and 0.9531/?i, respectively. The other parameters 
are the same as in Fig. 2. The exact results of the propagation 
constants for the composite waveguides are also included in 
Fig. 5 for comparison. One sees that the results obtained from 
Eq. (21) are in good agreement with the exact results. 

Finally we analyze qualitatively the total energy flow of 



Fig. 5. Propagation constants f}± versus the distance d between the core 
layers for the composite waveguides with various NIM layer thicknesses. The 
NIM waveguide in isolation guides a single mode forward (P22 > 0) for all 
cases. The lines represent the results obtained from coupled mode theory, and 
the marks represent the results calculated from Eq. (22). 



the modes in the composite waveguide, which is given by 
S = U + PU, where the superscript " + " denotes the con- 
jugated transpose of matrix. For the waveguide system with 
P22 > 0, P is a positive Hermitian matrix, which leads to 
S > 0, indicating that the modes of the composite waveguide 
are always forward modes, i.e., the total energy flows of the 
modes are codirectional with their phase flows. Therefore, the 
guided modes shown in Fig. 5 are all forward modes. This 
is verified by the exact calculation of the total energy flows 
of the modes in the composite waveguide. On the other hand, 
for the waveguide system with P22 < 0, P is an indefinite 
Hermitian matrix, and it can be written as P = R + JR, where 
J = diag[l,-l], and R+ = [ri,r 2 ] diagflyW vM, 
with A„ and r n (n = 1, 2) being the eigenvalues and 
orthogonal eigenvectors of P, respectively. Let V = RU = 
col[«i, V2], then the coupled mode equations are expressed as 
dV/dz — jTV, where T is an anti-Hermitian matrix, given 
by T = J[R + ]~ 1 (PB + KJR^ 1 = J[R + ]~ 1 HR~ 1 , with 



T — 

(\V1/V2? 

find 



T 2 \. We also have S = V+JV 



V2 



1 ) I 1 • Solving the coupled mode equations, we 



P± = 



{Tn + T 22 ) ± y/{Tu -T22) 2 -4|Ti 2 | 



and 



\vl|v2\ = \t±^fe^\\ 



(23) 



(24) 



where £ = | (T tl - T 22 ) /2T 2 i| > 0. Note that it is not 
necessary to take the same sign simultaneously in Eqs. (23) 
and (24). As seen from Eq. (23), when £ < 1, the propagation 
constants of two modes in the composite waveguide have an 
imaginary part. In this case we have |t>i/v 2 | = 1 and then 
5 = 0, which indicates that the modes in the composite 
waveguide are indeed evanescent waves, though they have a 
nonzero real part of the propagation constants. In the case 
when £ > 1, /3± are both real, i.e., the modes in the composite 
waveguide are both propagating waves. In this case one has 
\vi/v 2 \ > 1 and then S > 0, if we take " +" in Eq. (24), while 
\vi/v 2 \ < 1 and then S < 0, if choosing " - " in Eq. (24). 
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Thus, one of the two modes in the composite waveguide is a 
forward mode, and the other must be a backward mode, whose 
total energy flow is contradirectional with its phase flow. Since 
I-P12I, I-P21I *C |Pn|, P22 1 m general cases (i.e., except in the 
case of super strong coupling), the matrix P « diag[Pn, P22], 
thus we have approximately R « diag[^/|Pn|, y^jp^T]. 
Substituting this equation into the expression of T, we find 
that Tn »A. T22 ~ fa. So we infer that S > for (3+ and 
S < for 0_, if /3i > fo, otherwise S < for + and S > 
for /3_, if /3i < /?2- Therefore, in Fig. 4, the guided modes 
with propagation constant [3 + are forward modes, while the 
ones with propagation constant /3_ are backward modes. This 
is also verified by the exact numerical calculation of the total 
energy flow. 

V. Discussion and Conclusions 

The wave propagation and coupling in parallel planar NIM 
and PIM waveguides have been studied theoretically. The 
coupled mode equations for such a waveguide system has 
been developed, which has been verified through calculation 
of the propagation constants of the modes in the composite 
waveguide. It has been shown that if the NIM waveguide 
in isolation guides its mode forward, the properties of the 
NIM and PIM waveguide system are similar to those for 
a conventional waveguide system. However, in the case if 
the NIM waveguide guides its mode backward, interesting 
phenomena then appear. When the propagation constants of 
the individual waveguides are nearly equal, there exist only 
evanescent modes in the composite waveguide. The solution of 
the coupled mode equations for this case varies exponentially 
with the coupling length, which never happens for a conven- 
tional waveguide system. A coupler operating in this case is 
very easy to fabricate, as its coupling efficiency is insensitive 
to the coupling length. When the propagation constants (3% and 
P2 differs significantly, the modes in the composite waveguide 
may change from evanescent waves to propagating ones as d 
increases. In the latter case the solution of the coupled mode 
equations is also a periodic function of the coupling length, 
but the coupled power is output in a direction opposite to that 
of input power. Correspondingly, the modes in the composite 
waveguide are propagating waves, and one of the modes is a 
forward mode while the other is a backward one. Finally we 
should indicate that as the NIM medium is strongly dispersive, 
the various phenomena mentioned above may happen in the 
same waveguide system at different frequencies. 
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